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We define an effective action for spin foam models of quantum gravity by adapting 
the background field method from quantum field theory. We show that the Regge 
action is the leading term in the semi-classical expansion of the spin foam effective 
action if the vertex amplitude has the large-spin asymptotics which is proportional to 
an exponential function of the vertex Regge action. In the case of the known three- 
dimensional and four-dimensional spin foam models this amounts to modifying the 
vertex amplitude such that the exponential asymptotics is obtained. In particular, 
we show that the ELPR/FK model vertex amplitude can be modified such that the 
new model is finite and has the Einstein-Hilbert action as its classical limit. We 
also calculate the first-order and some of the second-order quantum corrections in 
the semi-classical expansion of the effective action. 



1. Introduction 



The problem of determining the classical limit of a spin foam model and the corresponding 
quantum corrections is one of the least understood and it represents the greatest obstacle 
for formulating a realistic spin foam model of quantum gravity. Certain clues have been 
obtained over the years, and the first one was the result obtained in [1] about the speed of 
convergence of the Euclidean Barret-Crane (BC) spin foam model state sum as a function 
of the edge amplitude. It was observed that a typical convergent BC state sum has the 
dominance of the small spin configurations. However, one can not obtain a sufficient 
information about the classical limit of a spin foam model from the partition function, 
and in [2] a study of the large-distance asymptotic behaviour of the graviton propagator 
for the BC model was started. It was discovered that the graviton propagator will have 
the correct large-distance asymptotics if the boundary spin-network wavefunction has a 

^Member of the mathematical physics group at the University of Lisbon. 



1 



certain Gaussian form. However, in [3] it was pointed out that the tensorial structure 
of the BC model graviton propagator is not correct. This problem was caused by the 
absence of the intertwiners in the BC model. Without intertwiners, one cannot construct 
the complete Hilbert space of Loop Quantum Gravity (LQG) on a three-dimensional 
boundary. 

Furthermore, it was pointed out in [4] that it is difficult to construct a physical spin- 
network wavefunction which has the Gaussian form used in [2] and it is not clear whether 
such a wavefunction exists. However, it was noticed in [4] that if a physical wavefunction 
has the Rovelli Gaussian form in the limit of large spins, then such an asymptotics is 
sufficient to ensure the correct large-distance asymptotics for the graviton propagator. 
Further studies on the large-spin asymptotics of a class of physical wavefunctions which 
can be constructed as three-dimensional spin-foam state sums have revealed that their 
large-spin asymptotic form can be a Gaussian function, but it can not be of the Rovelli 
type, sec [5, 6]. The source of this problem was identified in the large-spin asymptotics of 
the ELPR/FK spin-foam vertex amplitude [6]. Namely, this asymptotics is of the form 

n) ~ — , (1) 

where j are the spins of the faces meeting in a given vertex, n are the corresponding 
coherent state vectors and SyR{j,n) is the corresponding 4-simplex Regge action, see 
[8, 9]. V{j) is a homogeneous function of order 12, N± are homogeneous functions of j 
of order zero and a is a constant. On the other hand, it was shown in [6] that a vertex 
amplitude whose asymptotics would be given by (1) such that = or A^_|_ = would 
give the correct graviton propagator asymptotics. 

This result can be easily understood from the path-integral point of view, since the 
vertex asymptotics with N_ — gives the state sum which for large spins looks like the 
usual path-integral 

Z ^ jV(j) e*^^'^) , (2) 

where S{(j)) is the Einstein- Hilbert action. On the other hand, the asymptotics (1) with 
N+N_ ^ gives the state sum which for large spins looks like 



(3) 



and this is an unusual form of the path integral. 

The authors of [10] argued that the path integral (3) could still give the correct classical 
limit, if an appropriate boundary state is used. However, as argued in [6], it is difficult 
to prove that such a state exists. In this paper we will demonstrate that the usual path 
integral (2) gives the correct classical limit, and this will be done by using the effective 
action. At the heuristic level this can be seen from the background field method definition 
of the effective action 

^m/n ^ jj^^^is(^+m^ (4) 

see [11] for details. The definition (4) implies that S{(f)) is the classical hmit of r(0) in 
the case of the usual path integral, while in the case of the path integral (3) we obtain 
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It is clear that this expression will not give 5(0) as the classical limit of r(0). However, 
in order to have the state sum which for large spins looks like (2), we need to change the 
vertex amplitude such that the large-spin asymptotics is given by 



W{j,n) 



V{j) 



(5) 



In this paper we will show that this can be done in the case of the EPRL/FK class 
of spin foam models, and that the length-Rcggc action can be obtained in the large- 
spin limit. In section 2 we discuss the known approaches to the problem of the classical 
limit of spin foam models and conclude that the effective action approach based on the 
background field method is the most feasible. We give the definition of the effective 
action for spin foam models and in section 3 we study the case of the Ponzano-Rcggc 
model. We show that the corresponding effective action has the desired classical limit 
if the vertex amplitude is redefined such that the large-spin asymptotics is proportional 
to the exponential function of the vertex Regge action times the imaginary unit. We 
calculate the corresponding effective action in the large-spin limit and show that the 
Regge action is the dominant term in the semi-classical expansion. We also formulate a 
finite Lorentzian Ponzano-Regge model such that its classical limit is the Einstein-Hilbert 
action. In section 4 we study the case of the Lorentzian ELPR/FK spin foam model and 
show that the vertex amplitude can be redefined such that its large-spin asymptotics is 
proportional to the exponential function of a constant times the vertex Regge action. The 
corresponding effective action can be calculated in the large-spin limit by using the same 
techniques as in the three-dimensional case, and the result is the length-Regge action plus 
the quantum corrections. In section 5 we present our conclusions. 

2. Classical limit of spin foam models 

Spin foam models of quantum gravity in D = 3 and D — 4 spacetime dimensions are 
described by a partition function of the form 



where j = (ji, • • • , J/, • " " > Jf) and l — {li, ■ ■ ■ , li, ■ ■ ■ , ll) are labels of a spin foam whose 
two-complex a is dual to the simplicial complex obtained by triangulating the spacetime 
manifold. The faces / and the links / of a carry the labels jf and li, respectively, and 
these labels are irreducible representations and the corresponding intertwiners of the group 
Spin{D). Since Z is a complex number, it is not possible to extract the classical limit 
from it and one needs to analyze the boundary wavefunctions or to analyze the effective 
action. 

A boundary wavefunction ^(s) is associated to a boundary spin network s = {'y,jb, Lh)i 
where 7 is the boundary one-complex of a and (j, l) are the corresponding labels, ^'(s) is 
constructed from (6) such that the summation is restricted to spin foams whose boundary 
spin network is s. Therefore 



Z 



(6) 



3,1- f I 



(7) 



3,1- f I 
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where the amphtudes W are the same as the amphtudes W for the faces and the hnks 
not belonging to the boundary, while for the boundary faces and links there is a choice 
which ensures good gluing properties, see [12]. 

Note that the construction (7) gives just one boundary state 

\^) = T.^{s)\s), (8) 

s 

while we know from the canonical LQG that there are many different physical states. 
Especially important physical states are those which describe fiat or constant curvature 
spatial manifolds. Therefore the definition (7) has to be changed such that the information 
about the background triads Eo{x) is included, where a; is a spatial coordinate. In the case 
of Euclidean canonical LQG one can show that such a wavefunction has a form similar 
to (7), but the spin network s is replaced by a spin network s where s is s with edge 
insertions fii{Eo{l)), where £"0(0 — Ii Eo{x)dx, see [13]. The insertion functions can be 
chosen freely, and an appropriate choice arc the Gaussians centered around Eq{1). In the 
case of a flat geometry, all Eq{1) can be taken to be approximately the same, and the 
corresponding area of a triangle is proportional to Jq. This parameter will set the length 
scale, so that one introduces the insertions into the boundary spin network of (7) which 
will be functions of jo- 

Given such a ^(s,jo)i there is the corresponding connection wavefunction ^o(^)) 
which can be obtained by the loop transform. By writing 

'^q{A)^R{A) e^^(^)/\ 

one would have to show that 

s{A) = So{A) + nSi{A) + o{n^) 

where Sq{A) satisfies the Hamilton- Jacobi equation for canonical general relativity (GR). 
It is obvious that this is an extremely difficult way to obtain the classical limit. 

An easier approach would be to calculate the graviton correlation functions for the 
boundary state (8) with insertions, so that 

Gn{Xi, ■■■,Xn)^Yl %{s)^o{s'){s\h{xi) ■ ■ ■ h{Xn)\s') , 
s,s' 

where h is the graviton operator. This was the approach started by Rovelli [2], and it can 
be shown that G2 has the correct large-distance asymptotics if 

^o(s) « iVexp I -- J2 (^ii'Ui - Jo)iJi' - Jo) ' (9) 

for large spins, where a is a constant matrix [4]. 

However, the correct asymptotics of G2 does not guarantee that the classical limit 
of a spin foam model is general relativity. Namely, one has to show that all correlation 
functions Gn correspond to the ones for the EH action in the classical limit. This is 
equivalent to demonstrating that the effective action 

r[^] = H C„ / dXidti ■■■ dXndtnGn{xi,ti, ■■■,Xn, tn)h{xi, ti) ■ ■ ■ t^) , (10) 

n>2 '' ■' 
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has the classical limit which is given by the Einstein-Hilbert action, where {x, t) is a 
spacetime point coordinate and Gn is the extension of G„ when tk ti. The correlation 
function approach is less difficult then the wavefunction approach, but it still requires a 
lot of work. 

A simpler method to compute the effective action would be a method where F is given 
as a functional of the spacetime metric g rather then the functional oi h = g — rj, where rj 
is the fiat metric. The background field method (BFM) of computing the effective action 
in quantum field theory (QFT) [11], is a method convenient for this purpose. In the case 
of quantum gravity, the BFM approach suggests the following relation 

^ngyn ^ J^j^^.sig+h)/n^ (11) 

where S{g) is the Einstein-Hilbert action. The expression (11) is formal and has to be 
defined, and in the perturbative QFT approach it amounts to gauge fixing of the diffeo- 
morphism gauge symmetry and implementing a regularization/renormalization procedure. 
This has to be done because the theory allows arbitrary short distances and hence the 
infinities appear. Since GR is a non-renormalizible theory, the corresponding effective 
action can not be determined uniquely. 

In the case of spin foam models, the problems with QFT infinities are avoided because 
the theory has a natural short distance cut-off. Namely, the basic degrees of freedom 
are the SU{2) spins jf G N/2, and these are essentially the areas of the corresponding 

triangles, since Af (x Ip^Jjfijj + 1/2). Hence there is a short-distance cut-off of order of 
the Planck length Ip. There is no need for a gauge- fixing procedure, since jf are triangle 
areas, and these are diffeomorphism invariant. The only infinites which can appear are in 
the large jj region, which correspond to large- distance infinities (also known as the infra- 
red infinities in QFT), but these can be easilly dealt with, by introducing the appropriate 
negative powers of jf in the spin-foam amplitude, see [7]. 

The path integral (11) takes the following form in the case of spin foam models 

^ru,) ^j2Y[ W,{jf + j'f) n W,{ii + i\) n A(j^(,) + , + .;(,)) , (12) 

where the spin-foam two complex is closed, (j, i) is a configuration representing the back- 
ground or the classical values of the spin-foam labels, while the summation is over the 
fiuctuations (j', t') around the classical background (j, C). The calculation of F will sim- 
plify if we take that the background spins are large and that (j, i) is a stationary point of 
5'(j, i), where 

e^o-.o =nw^2(j/)nw^i(^on^o/M>^K^)) (13) 

is the partition function spin-foam amplitude. The latter condition is also used in the 
QFT version of the BFM approach, where it takes the form of background metric being 
a solution of the Einstein equations. 

As discussed in the introduction, the vertex amplitude A should be a function of the 
Wo amphtude such that A has the asymptotics 



for j ^ oo, where 

SvR=j:jf(^f (15) 

is the vertex Regge action and 9f are the dihedral angles, while Vp{j) should be a homo- 
geneous function of order p > 0. The role of the function Vp{j) is to ensure that the state 
sum (12) is convergent, see [7]. 

The requirement (14) is essential, since it will give the Regge action 

SR^j:jfSf (16) 

/ 

as the classical limit of the effective action. Namely, if we take the background spins to 
be large, then we can use the asymptotic formula (14) in the state sum (12), which then 
produces the Regge action in the exponent due to the identity 

Sn^Y.SvR + 27:Y,kfjf, (17) 
^ / 

where kf are integers. 

It will be important to notice that the deficit angle 5/ is given by 

Sf^27r-Y: O'f. , (18) 

vCf 

for a spacelike face /, where 9'^^ — n — 9fy is the interior dihedral angle for the 4-simplex 
ay dual to A spacclike / means that the dual triangle Ay is timelike and belongs to 
a^. When / is timelike, which means that A/ is spacelike, then 

Sf^E Of. , (19) 

VCf 

where 9fy — Qfy and ©/^ is the boost parameter between the normal vectors of the two 
tetrahedrons of a. which share the triangle A/, see [14]. 

4. The three-dimensional case 

Let us explore first the effective action given by (12) in the simpler case of three-dimen- 
sional (3d) spin foam models. Consider the Ponzano-Regge model partition function 

ZPR = E n(-l)''^dimiy n WUm) (-l)^-^(^)+-+^«(^) , (20) 
j f 

where W is the 6j-symbol, see [15]. The immediate problem with ZpR is that it is 
divergent, so that it has to be regularized. This can be done by introducing a maximum 
spin, or by dividing W by an appropriate power of the product of the dimensions of the 
vertex spins. The later regularization will be more convenient for our purposes. Either 
way one loses the triangulation independence, but we will see that this is not going to be 
a problem for our purposes. 
The next problem is that 
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for large spins, where V{j) is the volume of the vertex tetrahedron. According to our 
approach we then need to change the vertex amplitude such that the new asymptotics is 
given by (14). In order to achieve this consider 



w = Ww + Vvw^ - 1 . 

It is easy to show that W ~ e*'^"^ for large spins since (21) implies 



(21) 



W 



2VV 



Let us introduce a modified vertex amplitude 



W{j) 



(22) 



Then A{j) will have the asymptotic form (14) with p = 6p' + |, when all of the six spins j 
are large. The parameter p' has to be chosen such that the state sum (12) is finite. This 
can be done because W/y/V is a limited function. Then exists M > such that 



I^(J)I < 



M 



Consequently 



(23) 



(24) 



j f j f 

where Zp is the partition function associated to the state sum (12). The Uf denotes the 
number of vertices of a face / and since n/ > 2, we obtain the last inequality. The last 
sum in (24) will be convergent for p' > 1. Therefore Zp will be convergent for p' > 1. 
One can find a better estimate for the lower bound for p' by using a better estimate for 
W/VV, but for us it is important that such p' exist and that their values are independent 
from the triangulation. 

We arc interested in calculating the effective action when all the background spins in 
the state sum (12) become large. Then one can approximate each vertex amplitude in 
(12) by using (14), since all the spin labels j + j' are large. Consequently 



(25) 



j' f 



where we used ^2(7 + j') ~ 2(j + j') and pf = p' + 3/2. Note that the sign factors in 
the face and the vertex amplitudes, see (20), combine with the sum of the vertex Regge 
actions such that the Regge action is obtained in the exponent of (25). 

The main contribution in the state sum (25) comes from << jj, since the weights 
Uf +jfy''^f^f are maximal for = and drop-off as negative powers of j^. We can then 
use 



(j +/)-'" 



1-^ 



1 — m h mim + 1) — - + 

J 



which is valid for / < j. Consequently 



.iSR{.3+3') 



1 — m 



3'f 



3 f 



f — + mfimt + 1) .„ 
Jf Jf 



+ 



(26) 
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where m/ = n/p/ — lisa positive number. 

Let us choose the background spins jf such that they correspond to a stationary point 
of the Regge action SR{j). This is a standard procedure in the case of QFTs, and the 
idea is to simpUfy the calculation, since the stationary point restriction of the background 
spins does not affect the important features of the effective action. Therefore we can use 
the approximation 

SrU + f) ^ Snij) + IESr ffWfj'f , 
where S'^^ji^j) is the Hessian matrix for SR{j). Consequently 

e^m ^ Ne''nU)-J:fmfinj,j2e'is'^0)fjV2 Y[(l-mfi + --) , (27) 

j' /V jf J 

where {S'li{j)j'j') = J2fj'S'^ff,{j)j'fj'f,. The sum in (27) can be approximated by an 
integral over Xf — j'f/jf variables, and this integral will be given as a sum of the integrals 
of the following type 

/ d^X X^^ ■ ■ ■ Xp^ exp I - ^ Sji^^XjnXn I ■ 
\^ m,n ) 

These integrals can be calculated by taking the derivatives of the generating function 

l{ix) = / d^X exp ^ ^ S'^^^XmXn + ^ /^r, 



"mXm I ) 



m 



at /i = 0, where /(/x) is given by 



and ji^ = {fii, ■ ■ ■ , fip). 

This calculation can be simplified by using 

(1 + x)~^ — — Q-mx+mx^ /2+--- 

so that 



3' 

where 



Then by using (28) 



exp iiJ2f f ^f^f'~2'^~^ 
TV' exp I - ^ m/ log + iSnif) 1 ^ 1 5.. ^ (3°) 



/ 



(det5^^(j))V2 
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where G is the inverse matrix of S" . Consequently 

^ TV'exp [-Y^rriflogjf + iSr{j) - Irrlog^^^O') + ) ■ (31) 

V / 2 /,/' ^^f^i' j 

The equation (31) implies that 
m ^ SR{j)+iY.mf\og3f + l:TT\og~S'i,{j)-i\ogN' + Y.^fmf,^^ 

~ SrU) + iE^f log jf + '-Tr log 5^(i) - i log N' + 0(1/ j) , (32) 
/ ^ 

where 0{l/j) denotes the terms which scale as 1/j when jf — >■ jjf. When necessary, it 
will be understood that notation 0{j^) also includes the subleading terms. 

Note that Sr — 0{j), while log and Tr log 5"^ are of 0(log j). Therefore (32) implies 

r(Aj) ^ Aro(i) + (log A) ri(i) + v-,{j) + a-^ r3(i) + 0{i/\^) , 

in the limit A — ?> oo. The dominant term in the large-spin limit will be Fq, which is the 
Regge action (16). Hence we can say that the classical limit of the effective action is the 
Regge action, which is a discretization of the Einstein-Hilbert action. This means that if 
we start refining the spacetime simplicial complex, Tq will tend to the EH action. 

Note that the quantum correction terms in (32) arc imaginary numbers, while the 
effective action should be a real function. The same problem appears in QFT, since in 
that case 



The stationary phase approximation then implies 

r(0) ^ 5(0) +i-Tr log 5"(0) + 0{h^) 

for a stationary point 5"(0) = 0, so that T[(f)) is not a real action. This problem is resolved 
by resorting to the Wick rotation. Namely, by performing the Wick rotation t — ?> it, where 
t is the time coordinate, one passes to the theory in the Euclidean metric and defines the 
effective action which is real 



Consequently 



r(0) ^ 5(0) + I Tr log 5"(0) + 0{fi') . (33) 

Then in the Euclidean effective action (33) one replaces the Euclidean metric with the 
Minkowski metric, and obtains a real effective action. In the case of spin foam models, 
there are no spacetime coordinates and there is no background metric, so that one cannot 
perform the Wick rotation. However, note that the procedure used in QFT is a redefinition 
of a complex function F = To + iFi into a real function Fo + Fi. Therefore we can define 
a real effective action by 

r ^ ReT + ImT . (34) 
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The definition (34) is not unique, since one can also use Re T — Im T. This ambiguity can 
be only resolved by an experiment, but we will use (34), because it agrees with the QFT 
theory sign. Therefore for large spins we obtain 

m ~ Sr{j) + ^mf\ogjf + ^Tr log 5^(j) + 0(1/ j) . (35) 

We would like to make the following remarks. The partition function Zp, which corre- 
sponds to the modified vertex amplitude (22) is finite, but it is not a topological invariant. 
This is not a problem, since our goal is not constructing manifold invariants, but obtaining 
a quantum theory of gravity whose classical limit is general relativity. This is achieved 
by requiring that the state sum (12) is finite and that the classical limit of the effective 
action is the Rcgge action (16). Although our construction is triangulation dependent, it 
still leads to a topological theory in the continuum limit. Namely, if we start refining the 
triangulation, the Regge action will become the EH action, which defines a topological 
theory in three spacetime dimensions. 

The trace-log term in (35) is a discretization of the usual trace-log term from QFT. 
Since the QFT trace-log term is divergent, the spin foam version can be considered as a 
regularization of the QFT counterpart. In contrast, the mlogj terms in (35) do not have 
an analog in the QFT case, and their presence is a feature of the model. It is not clear 
what is the smooth limit of the mlogj terms and whether their presence is a good or a 
bad feature of the model, since we do not know experimentally what are the quantum 
gravity corrections. 

Note that one can define a spin foam model with a simpler vertex amplitude than (22) 

A(j) = ^ , . , . (36) 
nfe(dimjfe)f ^ ^ 

The corresponding state sum would look like the path integral for a Regge gravity model 
where the lengths of the edges of a triangulation are positive half-integers. One can enforce 
the triangle inequalities by inserting the dual-edge (triangle) amplitudes proportional to 
the theta spin network evaluation, so that 

^=En(-i)'^'^dimj^n^(Mo)n^o/M)(-ir^^^^---^^'^^''^- i^v 

j f IV 

The corresponding effective action will be also given by the expression (35). 

The spin foam model defined by the A amplitude can be easily extended to the 
Lorentzian case, simply by replacing the vertex Regge action S^r in (36) by its Lorentzian 
analog. We will label the edges of the triangulation with the unitary irreps jf from the 
discrete series of representations of Spin{l, 2) = SL{2, R), which makes sense if the edges 
are spacelikc. Then the deficit angle is given by the 3d analog of (19), so that the sum of 
the vertex Regge actions will be equal to the Regge action. Therefore we will not need the 
sign factors in the face and vertex amplitudes, which were necessary in the Euclidean case, 
in order to obtain the Regge action in the classical hmit of the effective action. Hence 
one can define a finite Lorentzian 3d quantum gravity spin foam model whose partition 
function is given by 

zl^y: n(2i/ + 1) n ^um) n MJm) . m 
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The corresponding effective action will given by (35) in the large-spin limit so that in the 
smooth spacetime limit one will obtain the EH action. 



5. The four-dimensional case 

We will consider the ELPR/FK spin foam models [16, 17], since this is the only class 
of four-dimensional (4d) spin foam models that has a well-defined LQG theory on a 3d 
boundary. The partition function is given by (6) such that W2{j) = dimj, Wi{j,i) = 1 
and 

ni>0,-,nB>Oo=l ''^ 

where 7 is the Barbero-Immirzi parameter, W15 is the 15j-symbol for the unitary repre- 
sentations (n, p) of the Lorentz group and / are the fusion coefficients, see [16] for the 
details. W2{j) was originally chosen to be a quadratic function [16], but it has been re- 
cently argued in [12] that the linear weight is more appropriate. In any case, the essential 
features of the effective action are the same. 
A more convenient form of Z is 

^ = E n^™.?/ n ^Uf{v),riif{v)) , 

where each Li in a spin foam from the sum (6) is replaced by four unit 3-vcctors nif. An 
n;/ is a 3- vector orthogonal to the triangle dual to face /, such that this triangle belongs 
to the tetrahedron dual to a link I. 

For the purposes of calculating the effective action wc only need to now the asymptotics 
of W{j, n) when all the spins j are large. This asymptotics is given by 

wu, n) ^ n^T^ ^^^^ 

where o; = 7 when the 4-simplex boundary geometry is Lorentzian, while a = 1 when 
the 4-simplex boundary geometry is Euchdean, see [9]. The real numbers N±[a) are 0(1) 
functions of the spins, while V — 0{j^'^). There are also degenerate configurations of 
spins for which W{j,n) f» N{j)/V{j) where N{j) — 0(1) and V{j) = 0(j^^), but their 
contribution to the state sum is negligible. Otherwise, W drops off faster than any power 
of 

In order to obtain a correct Fq we need to redefine the vertex amplitude. Let us 
introduce a new vertex amplitude W{j, n) such that 

W= (40) 

This formula follows from the following relation between W and W 

N^W + N_{W)-' 
V{j) 
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which ensures that 

for large spins. One can then define 

A(n n) = ^^-^'^^ (41) 

where p is sufficiently large such that Zp is convergent. Such a p, which does not depend 
on the triangulation, can be always arranged, as shown in [6]. 

The amplitude (41) has the desired asymptotics and it gives a finite partition function. 
We can now use the formula (12) for the effective action, which takes the form 

e^ro>i ^ ^ n[2(i/ + i}) + nUMJnv) + ■ (42) 

j',n' f V 

The labeling of the ELPR/FK model is consistent with a spacelike triangulation. In that 
case kf — in the formula (17), so that when (j'l, • • • , jV) (oo, • • • , oo) we obtain 

~ ^ H + j})i-f/'"/e^^«(^'+^"'"+"') , (43) 

where Pf > p and m/ is the number of vertices belonging to a face /. Note that the 
contribution to e^^ of the configurations {j+f, n+n') which are not geometric is negligible 
compared to (43), since A{j + j',n + n') decreases exponentially with large spins. The 
contribution of degenerate geometric configurations is also negligible, since one sums over 
a lower-dimensional sub-space in the space of spins. 

Let j and n be the background spin foam labels such that (j, n) is a stationary point 
of Sji{j,n) and all the 4-complexes have the Lorentzian geometry. Then we can use the 
formulas from the 3d case, and we obtain 



e 



imn) ^ AT ^ g-E/^/l°s(i/+i/0+iSfl(j,n)+i(S^^^i'j'+2S^^^„i'n'+S^„„n'n'> 
j',n' 



j',n' 

where Cf = PffUf — 1 and 

By performing the Gaussian integrations in (44) we obtain 

gir(i,n) ^ ^i^-Y,fCf\og{jf)+iSR{j,n)-\Tr\ogS'^{j,n)+\Yjfj'^f^f'(^RSf'ii^^) ^ ^45^ 

where 

'-'R I q/f q/i j 1 

\ '-'Rjn '-'Rnn / 

Gr fft is an element of the jj block of the matrix {S'^)~^ and Cf = Cf/jf. 
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The equation (45) implies 
m n) « Sn{j. n) + Y.^f log Jf + \Tr log n) +Y. CfCf, ^MlMl , (46) 

where we have used T ^ ReV + ImV and we have omitted the constant logA^'. As in 
the 3d case, the dominant term is Sr, which is of while the other terms are of 

subleading orders, namely of O(logj) and 0(l/j), respectively. 

The equation (46) implies that the classical limit of the ELPR/FK effective action is 
the area-Regge action 

SR{3,n)=^Y.3A3f,n). (47) 
/ 

However, we can require that the background spin foam (j', n) is also a stationary point of 
the partition function amplitude, sec (13). This requirement imposes further restrictions 
on the background spin foam, and it was shown in [18, 19] that such stationary points 
correspond to Regge geometries when jf are all large. This means that there is an as- 
sigment of the edge lengths L^, satisfying the triangle inequalities, such that jf oc Af{L) 
and n — n{L), where ^/ is the area of the triangle dual to a face /. 

More precisely, the stationary point equations imply that the four normals n associated 
to a tetrahedron r, which are determined by the four-geometry of a four-simplex a which 
contains r, are the same as the four normals n' determined by the four-geometry of another 
a' containing r. Consequently, the length of an edge e of r is the same when calculated 
from the ten triangle areas jf of cr or when calculated from the ten triangle areas j'f of 
cr', so that 

L^:\3i,---,3i.)^L^f\3[,---,j[,)- (48) 

This is precisely the constraint which turns an area-Regge action into a length-Regge 
action, see [20]. Hence 

SR{3.n) = lY.JfWf{L) ^^,Y.Mm{L) , (49) 
/ / 

where the labels L denote the lengths of the edges of the triangulation and Ip is the Planck 
length. 

If we refine infinitely the spacetime triangulation, then the Regge action in (49) will 
become the Einstein-Hilbert action S'^j^/lGTr/p. Therefore the ELPR/FK spin foam model 
with a modified amplitude (41) will give general relativity in the limit of large spins and 
smooth spacetime. 

Note that the structure of the quantum corrections in the effective action (46) is the 
same as in the 3d case. 

6. Conclusions 

We have demonstrated that it is possible to construct spin foam models of quantum gravity 
in 3 and 4 spacetime dimensions such that the corresponding effective action has general 
relativity as its classical limit. The effective action method resolves the long-standing 
problem of how to compute the semiclassical limit of a spin foam model. The spin foam 
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model defined by the vertex amplitude (41) is the first example of a finite four- dimensional 
spin foam model with the correct classical limit. 

In 3d case we have shown how to regularize the PR model and how to modify the 
vertex amplitude in order to obtain the correct classical limit. The same can be done in 
the case of the Turaev-Viro (TV) model [21], which is a quantum group regularization 
of the PR model. The TV model effective action can be calculated since the large-spin 
asymptotics of the quantum 6j-symbol is known [22] . Because the vertex asymptotics is 
also of the cosine type, this means that the TV model vertex amphtude has to be modified 
in order for the effective action to yield the EH action with a cosmological constant in the 
classical limit. The same applies to the quantum group regularization of the ELPR/FK 
model [23, 24], since the large spin asymptotics will be a deformation of the classical group 
asymptotics (39). We were also able to construct a Lorentzian version of the PR model 
such that the corresponding effective action has the desired classical limit, see (38). 

The structure of the quantum gravity corrections in 3 and 4 dimensions is the same; 
however, their spacetime interpretation is different. In 3d we have jf oc so that the 
trace-log term in the 3d effective action is the discretization of the quantum field theory 
trace- log term, while the log j/ terms do not have the QFT analog. In 4d case, jf oc Af{L^) 
so that the hessian S'^j^, is not the same as the discretized version of the QFT hessian, 
which is S^ee'- Consequently the 4d effective action will have additional terms of the type 
log f{Lc/lp) where / is a homogenous function of order one. It would be interesting to 
analyze the implications of the log[Af{L^)/lp] and log f{L^/lp) terms for cosmology. 

Note that the terms proportional to the Gff matrix in (32) and (46) constitute second- 
order quantum corrections, which are of 0{l/j). In order to obtain all second-order terms 
we need to know the 0{l/j) correction terms to the vertex Regge action in the vertex 
amplitude asymptotics. 

It is not difficult to see that our formalism gives the area- Regge action as the classical 
limit of the appropriately modified Barret-Crane (BC) model, because the length con- 
straint (48) is not enforced. Hence this justifies the concerns raised about the BC model, 
which were based on the convergence speed of the partition function, see [1], as well as 
the concerns related to the graviton propagator for the BC model [3]. 

Note that our method can be applied to any state sum model of quantum gravity, 
which means that one can study the effective action for the 4d Regge model. 
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